Seeing the light : experimental signatures of emergent electromagnetism in a quantum spin ice 
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The "spin ice" state found in the rare earth pyrochlore magnets Ho2Ti207 and Dy2Ti207 offers a beautiful 
realisation of classical magnetostatics, complete with magnetic monopole excitations. It has been suggested 
that in "quantum spin ice" materials, quantum-mechanical tunnelling between different ice configurations could 
convert the magnetostatics of spin ice into a quantum spin liquid which realises a fully dynamical, lattice- 
analogue of quantum electromagnetism. Here we explore how such a state might manifest itself in experiment, 
within the minimal microscopic model of a such a quantum spin ice. We develop a lattice field theory for this 
model, and use this to make explicit predictions for the dynamical structure factor which would be observed in 
neutron scattering experiments on a quantum spin ice. We find that "pinch points", which are the signal feature 
of a classical spin ice, fade away as a quantum ice is cooled to its zero-temperature ground state. We also make 
explicit predictions for the ghostly, linearly dispersing magnetic excitations which are the "photons" of this 
emergent electromagnetism. The predictions of this field theory are shown to be in quantitative agreement with 
Quantum Monte Carlo simulations at zero temperature. 

PACS numbers: 75.10.Jm 75.10.Kt, 11.15.Ha, 



I. INTRODUCTION 

The idea that a strongly interacting quantum magnet might 
support a spin liquid phase which remains disordered even at 
zero-temperature has fascinated — and frustrated — physi- 
cists ever since the seminal "resonating valence bond" (RVB) 
paper of Anderson in 1973' . Such a phase, it was argued, 
need not support the spin waves found in conventional mag- 
nets, but could instead exhibit "spinons" with fractional quan- 
tum numbers. Forty years later, the search for quantum spin 
liquids goes on, but with strong grounds for encouragement : 
a growing number of quantum magnets have been identified 
which do not order down to the lowest temperatures mea- 
sured, many of which have low-temperature properties which 
hint at spinons-"-. At the same time, the "spin ice" materi- 
als Ho2Ti207 and Dy2Ti207 have emerged as text-book ex- 
amples of classical (i.e. entropy-driven) spin liquids^"— These 
highly-frustrated magnetic insulators show algebraic corre- 
lations of spins over macroscopic distance&2ri£ and support 
magnetic monopole excitations which provide classical ana- 
logues to the spinons envisaged by Andersor^ii^. 

Recently, the idea of a "quantum spin ice" has also attracted 
considerable interest. The family of rare earth pyrochlore s to 
which Ho2Ti207 and Dy2Ti207 belong includes other sys- 
tems in which quantum effects play a much more important 
role^. Perhaps the most widely studied system of this type 
is Tb2Ti207. Like the classical spin ices, the magnetism of 
Tb2Ti207 is controlled by the competition between strong 
Ising anisotropy, and dipolar interactions which are ferromag- 
netic on nearest-neighbour bonds, so it is expected to be an 



"ice". However, in Tb2Ti207, anisotropic exchange interac- 
tions also play an important role, and endow the spins with 
dynamicsi^^. A diffuse, liquid-like structure is observed in 
neutron scattering for a wide range of temperatures, with no 
conventional magnetic order observed down to 50mK, despite 
the fact that the typical scale of interactions between spins is 
closer to llK^^a^. Muon spin rotation experiments, mean- 
while, suggest that spins continue to fluctuate down to the 
lowest temperatures*^, and the most recent quasi-elastic neu- 
tron scattering experiments find evidence of power-law spin 
correlations at SOmK^^. Taken together, these facts make 
Tb2Ti207 a prime example of a three-dimensional, quantum 
spin liquid. 

The magnetism of Yb2Ti207 has also proved very inter- 
esting, with neutron scattering finding no evidence of order 
at temperatures above 210mK, and evidence for frustrated, 
anisotropic exchange interactions favouring significant dy- 
namics within an "ice-like" manifold of states Compa- 
rable studies of Pr2Sn207 suggest that it does not order down 
to 500mK, but with spins continuing to fluctuate^"—. There 
is also reason to believe that other Pr metal oxides, including 
Pr2Zr207, may prove a worthwhile hunting ground for quan- 
tum spin liquids^—. And, while the dynamics of the "clas- 
sical" spin ices Ho2Ti207 and Dy2Ti207 become very slow 
at low temperatures, neither system has ever been observed 
to order, despite the fact that the dipolar interactions present 
in these systems are expected to favour an ordered state 
All of this begs the question of how the classical spin liquid 
found in spin ice might evolve into a quantum spin liquid as 
quantum effects become more important ? 
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FIG. 1: (Color online). Spin correlations in a spin ice, as measured 
by quasi-elastic neutron scattering : (a) Correlations within the clas- 
sical spin ice configurations, showing the characteristic "pinch point" 
singularities, (b) Correlations in a quantum ice at T = 0, showing the 
suppression of pinch points by quantum fluctuations, (c) Correlations 
in a quantum ice at an intermediate temperature T = ca^^, showing 
how pinch points are progressively restored by the thermal excitation 
of magnetic photons. In all cases, results are shown for equal-time 
structure factors in the {h, h, 0) plane, for a polarised neutron scatter- 
ing experiment in the spin-flip channel considered by Fennell et al '^. 
Temperature is measured in units where c is the speed of light associ- 
ated with magnetic "photon" excitations, ao the lattice constant, and 

h = kn ~ 1- 



In fact spin ice is just one example of a much broader class 
of systems which obey the "ice rules". First introduced by 
Bemal and Fowler in 1933 to describe the correlations of pro- 
tons in water ice^^*^, the ice rules have since found application 
in models of frustrated charge orderi^Mi, proton bonded fer- 
roelectrics^ and dense polymer melts^. All of these systems 
possess a local "two-in, two-out" constraint, which can most 
conveniently be written in terms of a zero-divergence condi- 
tion on a notional magnetic field 



V -B = 0. 



(1) 



In the case of spin ice, B has the physical meaning of the 
local magnetisation of the system, and we can associate a 
field Bi with each spin on the lattice. For this reason, spin 
ice offers a beautiful realisation of classical magnetostatics, 
with local violations of the ice rules entering as point mag- 
netic charges (magnetic monopolesii^) and spin correlations 
which exhibit "pinch point" singularities in k-space 



(S'^(-k)S'^(k))classical OC [S^ 



(2) 



[Fig.fTJa)] corresponding to algebraic (dipolar) correlations in 
real spaca^ ''"'"^' . Since the ice rules can be satisfied by an ex- 
ponentially large number of proton (spin, charge, polymer. . . ) 
configurations^, they explain the residual entropy observed 
in both water ice^ and spin ice^ at low temperatures. Given 
this enormous reservoir of entropy, both spin ice and water ice 
are natural places to look for a quantum liquid ground state. 

The key ingredient needed to convert a classical ice into 
a quantum liquid is tunnelling between different ice config- 
urations [Fig. 111. This opens the door to a "quantum ice" : 
a unique, quantum mechanical ground state, formed through 
the coherent superposition of an exponentially large number 
of classical ice configurations. Such a state could have a van- 
ishing entropy at zero temperature, and so satisfy the third law 
of thermodynamics, without sacrificing the algebraic correla- 
tions and fractional excitations (magnetic monopoles) associ- 
ated with the degeneracy of the ice states. If reahsed in a spin 
ice, it would provide a concrete, three-dimensional example 
of the long-sought quantum spin liquid. 

Precisely this scenario was proposed by Moessner and 
Sondhi in the context of three-dimensional quantum dimer 
models^, by Hermele, Balents and Fisher in a quantum, ice- 
type model derived from an easy-axis antiferromagnet on a 
pyrochlore lattice^, and by Castro-Neto, Pujol and Fradkin 
in the context of a simplified model of water ice^. All of 
these models included tunnelling between ice (or dimer) con- 
figurations of the type illustrated in Fig. |2l In a spin ice, the 
dominant tunnelUng process involves flipping loops of spins 
which point nose-to-tail on an hexagonal plaquette, and the 
resulting dynamics are described symbollically by 

^tunnelling = -ff ^ [| O) (O | + | O) (O |] (3) 

o 

where g is the strength of the tunnelling matrix element, and 
^tunnelling acts on the space of all possible ice (or dimer) con- 
figurations. 
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FIG. 2: (Color online). An illustration of the simplest tunnelling 
process between different spin-ice configurations. The ice rules dic- 
tate that each tetrahedron within the lattice has two spins which point 
"in", and two which point "out". Where these spins form a closed 
loop on a hexagonal plaquette — here shaded red — the sense of 
each spin within the loop can be reversed to give a new configuration 
which also obeys the ice rules. 



Both Moessner and Sondhi^- and Hermele et alM also in- 
troduced an additional control parameter /i to the Hamiltonian 



tunneling 



6n, 



where 



6n, = mE[I(^>(«I + I^^>(^I] 



(4) 



(5) 



o 



This makes it possible to fine-tune the model to an exactly sol- 
uble Rokhsar-Kivelson (RK) point g = ^, where the ground 
state wave is an equally-weighted sum of all possible ice 
(dimer) configurations-?. The authors then argued, by con- 
tinuity, that a quantum liquid phase would occur for a finite 
range of parameters /i < 1 bordering on the RK point"*^ "*^ . 

The most striking feature of this quantum liquid is "light". 
Attempts to construct models with "artificial light" — gap- 
less photon excitations of an effective, low-energy U{1) gauge 
field — have a long history^. In recent years, it has been re- 
alised that large families of lattice models could, in principle, 
be described by such theories. These include abstract mod- 
els of "quantum oidef'^^^, Bose-Hubbard models bordering 
on superfluidity^, systems of screened dipoles^'*, and suitably 
adapted sigma models^. Reviews of these ideas can be found 
in papers by Montrunich and Senthil^ and Wen and Levin^. 

The way in which "light" arises in three-dimensional quan- 
tum ice and quantum dimer models is particularly simple. The 
ice-rules constraint Eq. ([T]) is most conveniently resolved as 



B{r) = V X A{r) . 



(6) 



The new feature which enters where there is tunnelling be- 
tween ice configurations is the fluctuation in time of the gauge 
field A{r). In conventional electromagnetism, this gives rise 
to an electric field 



£{r) 



dAjr) 



(7) 
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FIG. 3: (Color online). Ghostly magnetic "photon" excitation as 
it might appear in an inelastic neutron scattering experiment on a 
quantum spin ice realising a quantum ice ground state. The photon 
dispersion i^(k) is taken from lattice gauge theory developed in Sec- 
tion [nc] of this paper, convoluted with a Gaussian representing the 
finite energy resolution of the instrument. The intensity of scattering 
vanishes as / oc aj(k) at low energies. 



Then, following the heuristic arguments of Moessner and 
Sondhi"*- — or the more microscopic derivation of Hermele 
et alM- — it is reasonable to suppose that a quantum liquid 
found bordering the RK point ("H^ [Eq. Q] with /x < 1), 
would be governed by the Maxwell action 



1 

8^ 



^Maxwell — I dtSv 



£{rf - c^B{v) 



(8) 



Any state described by iSiviaxweii [Eq. (H))] automatically sup- 
ports linearly-dispersing transverse excitations of the gauge 
field .4(r) — "photons", with a speed of "light" c. On the lat- 
tice, such a magnetic photon would have a dispersion w(k) of 
the form illustrated in Fig. [3] 

Moreover, the fact that the spins now fluctuate in time, as 
well as space, introduces an additional power of k in energy- 
integrated (i.e. equal time) spin correlationsiid^, 



(S'p(-k)S'i.(k))quantum OC fc (5^ 



(9) 



which serves to eliminate the pinch points seen in quasi-elastic 
neutron scattering experiments [Fig. [lib)]''** . More formally, 
this theory is a compact, frustrated U{1) gauge theory on a di- 
amond lattice, and we will refer to the liquid state it describes 
as the quantum U{1) liquid below. 

The degree of fine-tuning in these arguments, and the need 
to introduce additional parameter ji [Eq. Q], might seem 
to render them of purely academic interest. However the 
idea of a quantum U{1) liquid found strong support in finite- 
temperature quantum Monte Carlo simulations of an ice-type 
model of frustrated charge order on the pyrochlore lattice^. 
Subsequently, it has proved possible to determine the ground 
state phase diagrams of both the quantum dimer model of 
Moessner and Sondhi^, and the quantum ice model of Her- 
mele et al^, from zero-temperature quantum Monte Carlo 
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FIG. 4: (Color online). Zero-temperature phase diagram of the 
model of tunnelling between ice states "H,! [Eq. [T6j, as determined 
by quantum Monte Carlo simulation in 58. The "quantum ice" point, 
fj. = 0, lies deep within a quantum liquid phase with low-energy exci- 
tations described by a lattice analogue of quantum electromagnetism 
This extends from a "squiggle" ordered phase, found for /i < —0.5g, 
to the exactly-soluble RK point i-i — g. (Here g is the strength of tun- 
nelling between ice states). 



simulation o ' • . Both models contains extended regions of 
a quantum liquid phase, connecting to the RK point. In both 
cases, this quantum liquid has low energy excitations which 
are described by a lattice analogue of quantum electromag- 
netismJMft^''. Significantly, in the case of the quantum ice 
model, this quantum liquid phase encompasses the "physical" 
point of the model /i = 0, and so does not require any fine- 
tuning [Fig.m^^. 

The theoretical possibility of a three-dimensional spin- 
liquid state with excitations described by a lattice analogue 
of quantum electromagnetism is now well-established. What 
remains is to connect these ideas with experiments. The pur- 
pose of this paper is therefore to set out predictions for the 
correlations which would be measured in neutron scattering 
experiments, ;/ such a state were realised in a spin-ice ma- 
terial. For concreteness, we work with the minimal lattice 
model introduced by Hermele et alM; transcribed to coordi- 
nates appropriate for a spin ice. More realistic generalisations 
of this model will be considered elsewhere^. 

In Section of the paper, we develop the mathematical 
formalism needed to describe the spin correlations and low- 
energy spin excitations in a spin ice with a quantum U{1)- 
liquid ground state. Using this theory, we make predictions 
for the photon dispersion CL>(k) and dynamical structure factor 
S'"^(k,aj), which would be measured in neutron scattering 
experiments. 

In Section Uni we make explicit comparison of the predic- 
tions of this theory with zero-temperature Quantum Monte 
Carlo simulations of the minimal, microscopic model of a 
quantum spin ice with tunnelling between different ice config- 
urations, Ti,^ [Eq. (fT6l)l. We find essentially perfect, quanti- 
tative agreement between simulation results and the field the- 
ory solved on a finite-size lattice, for a range of parameters 
< M < .9 which interpolate from the minimal model of a 
quantum spin ice (/^ = 0), to the classical correlations of the 
RK point ii^L = g). This analysis reinforces the conclusions 
reached in [58] about the existence of a quantum U (l)-Iiquid 
in this model, and puts the field-theory description on a quan- 
titative footing. 

In Section |IV] we make predictions for neutron scattering 
experiments carried out at finite temperature. In particular 



we analyse the way in which the characteristic "pinch point" 
structure in quasi-elastic scattering is lost as the system is 
cooled towards its zero-temperature ground state. We con- 
clude that the loss of the pinch points coincides with the pro- 
gressive loss of the Pauling entropy as the system cools into a 
unique, quantum coherent, liquid ground state. Thus the sig- 
nature features of the ice problem: pinch points and the Paul- 
ing entropy die together at low temperatures. We also give 
a brief discussion of the uniform magnetic susceptibility and 
heat capacity, in the low temperature quantum regime. 

Finally, in Section|V]we conclude with a discussion of some 
of the remaining issues relating to experiment. As far as pos- 
sible, each section of the paper is written so as to be self- 
contained. Readers uninterested in the mathematical develop- 
ment of the theory are therefore invited pass directly to Sec- 
tion|lIl]and Section|lV] referring to SectionlUas required. 



II. FROM QUANTUM ICE TO QUANTUM 
ELECTROMAGNETISM 

At first sight, an assembly of magnetic ions on a lattice does 
not look like a promising place to search for a gauge theory 
which perfectly mimics quantum electromagnetism. However 
in the simplest microscopic model for quantum mechanical 
tunnelling between spin configurations obeying the "two in, 
two out" ice rule, this is exactly what happens^i^^iS. In what 
follows we retrace the steps which lead from a spin ice system 
to a theory of electromagnetism on a lattice. 

In Section III Al we review the relevant microscopic mod- 
els. In Section HTbI we show how a lattice gauge theory re- 
sembling electromagnetism arises in these problems, recast- 
ing the earlier field-theoretical arguments of Hermele et al— 
in terms appropriate for a spin ice. In Section HTCl we explic- 
itly construct the magnetic "photon" excitations of this lattice 
gauge theory. In Section III PI we use the mapping between 
spins and photons to calculate the correlations between spins 
in a quantum spin liquid described by this lattice gauge the- 
ory. Throughout this analysis we set ?i = fc^ = 1, restoring 
dimensional factors only where we quote a result for the speed 
of Hght. 



A. Spins on a pyrochlore lattice 

The materials which we will seek to describe have mag- 
netic ions which a) have a crystal-field ground state which is 
a doublet, and b) occupy the sites of the pyrochlore lattice 
shown in Fig. |5l In the case of the spin ices Ho2Ti207 and 
Dy2Ti207, this doublet has Ising character (rare-earth mo- 
ments point into, or out of, the tetrahedra which make up 
the lattice), and the dominant interactions between these Ising 
spins are dipolarS. However, since these dipolar interactions 
are effectively self-screened, the correlations present in spin- 
ice are extremely well described by models with only nearest- 
neighbour interactions between spins^i^^i^i^i^. This approxi- 
mation gains further justification in "quantum spin ice" mate- 
rials such as Yb2Ti207, where magnetic moments are smaller 
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FIG. 5: (Color online). Structure of the pyrochlore lattice realised by the magnetic ions in spin-ice materials, a) The lattice is built of comer 
sharing tetrahedra, and can be decomposed into a set of A-sublattice tetrahedra (here coloured red) and B-sublattice tetrahedra (here coloured 
black), each of which forms an FCC lattice in its own right. The primitive unit of the pyrchlore lattice consists of a single tetrahedron with 4 
lattice sites. However it is also possible to define a cubic unit cell, of side ao, containing 16 lattice sites, b) Bipartite, diamond lattice formed 
by the centres of tetrahedra which make up the pyrochlore lattice. The bonds of this diamond lattice define the easy axes for spins in a spin 
ice, and play an important role in the lattice gauge theory of its excitations. 



than for Ho2Ti207 and Dy2Ti207, and exchange interactions 
play a much larger role. 

As a starting point, we can therefore consider the Hamilto- 
nian for a (pseudo) spin- 1/2 degree of freedom on a pyrochlore 
lattice, with the most general nearest-neighbour exchange in- 
teractions allowed by symmetry^ 

ns=i/2 = E {j^^-^m ^±(5+57 + SrS+) 

+J±± h,S+S++7*SrS7] 

+J,± [Sf (0,S+ + C* ST) + 1^]]] (10) 

Here we have followed the notation of Ross et al.—, in which 
the Sf is aligned with the local trigonal axes of the pyrochlore 
lattice on each site i, and 7^ and Qj are 4 x 4 complex uni- 
modular matrices encoding the rotations between these lo- 
cal coordinate frames. In the "quantum spin ice" Yb2Ti207, 
where the ground state doublet of Yb has XY-character— 
Eq. ( fTOb gives a good account of diffuse structure observed in 
neutron scattering experiments provided that the exchange in- 
teractions J±, Jz± and J±± are taken into account-i^i^. It also 
gives an excellent description of spin wave excitations about 
the saturated state of Yb2Ti207 in applied magnetic field, with 
parameters J^. = 0.17 ± 0.04 meV, J± = 0.05 ± 0.01 meV, 
J^± = -0.14 ± 0.01 meV and J±± = 0.05 ± 0.01 meV ob- 
tained from fits to data-^. The phase diagram associated with 
■Hs=i/2 [Eq. Cni] is explored in ll27ll35ll70ll . 

We can further simplify the problem by setting J±± = 0, 
and focusing on the limit Jzz >> J±,Jz± > 0. In this 
limit, the role of Jzz is to enforce the "ice rules" constraint, 
while J± generates dynamics, and Jz± lifts the degeneracy of 
ice-rule obeying states. Performing degenerate perturbation 
theory in the basis of (spin) ice configurations, and dropping 
terms which lead only to a constant energy shift, leads to the 



effective Hamiltoniansi 

"Heff = HtLnnelling + ^Ja (1 1) 

with 

■Htunnelling = ^5 X! [^fS^ S;;[S4 S^Sg" + h.c] (12) 
O 

where J2o '^^^^ hexagonal plaquettes in the py- 
orchlore lattice [cf. Fig|2l with 

12Ji 

.9=^, (13) 

and 

^J3 = -^3 E S^S| (14) 

where rui^^ over third-neighbours bonds (parallel to 

the nearest-neighbour bonds), with 

Js = ^ > (15) 

We note that, by construction, the Hamiltonian Eq. (fTTI) 
acts only on spin configurations satisfying the ice rules. This 
implies that, in performing the degenerate perturbation the- 
ory, virtual excitations of magnetic monopoles have been pro- 
jected out of the problem. This approximation will have little 
effect on the conclusions drawn in this paper, and could in 
principle be relaxed. 

It is also important to note that these spin ice configurations 
may possess a non-zero net magnetisation M. The tunnelling 
term 'Htunnelling [Eq. ( fT2b l generates dynamics by performing a 
cyclic exchange of spins on a hexagonal plaquette [cf. Fig.|2l. 
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This tunnelling process can be written symbolically as acting 
on a closed loop of spins [cf. Eq. Q]. Under these dynamics, 
the total magnetisation M is a conserved quantity. 

We make the final simplification of neglecting Hj^ 
[Eq. (fT4l )l and focusing exclusively on the spin-liquid 
favoured by the tunnelling term Htunneiiing — Eq. (fTZb or, sym- 
bolically, Eq. ([3]). The neglected term T^j, favours the six 
ice states with the maximum possible magnetisation per site 
m = (±l/-\/3, 0, 0) X S, etc., where S is the moment of the 
magnetic ion. We have confirmed through zero-temperature 
quantum Monte Carlo simulation of Heff [Eq. (fTTl il. that the 
system remains in quantum U{1) liquid ground state up to a 
value of J3 « 0.275, at which point it undergoes first-order 
transition into this ordered, ferromagnetic state. These results 
will be discussed elsewhere^. We note that a gauge mean- 
field theory for Hs=i/2 [Eq. (fTOb l predicts an intermediate 
"Coulombic ferromagnet" phase, in which the quantum C/(l) 
liquid spontaneously acquires a finite magnetisation for any 
finite J^jjaiiZS. This does not appear to be a ground state of 
the effective model "Heff [Eq. (fTTlil. 

Following Hermele et al^'' , it is useful to augment the min- 
imal model ^tunnelling with an additional, artificial, interaction 
term SHfj, [Eq. ©J. This renders the model exactly soluble for 
fj, = g. Thus the most general microscopic model we consider 
in this paper can be written symbolically as 



= -.9E0f^>(«l + |o>(o|] 

o 

o 

where "H^ acts on the space of all possible (spin) ice config- 
urations. This Hamiltonian is known to support a quantum 
U{1) liquid ground state for — 0.5.g < ^ < g^. 

It is important to note that this effective description of tun- 
nelling between ice configurations might equally have been 
derived for the model of hardcore bosons on the pyrochlore 
lattice considered by Banerjee et alM 



^charge-ice = -tJ2(^h+^]'' 



(17) 



with F > <. At 1/2-filling [{n) = 1/2], V selects 
charge configurations with exactly two bosons in each tetra- 
hedron of the lattice, and "Hcharge-ice is exactly equivalent 
to the pseudospin-1/2 model Eq. ( fTTI ). in the case where 
Jz± = J±± = 0. The leading tunnelling matrix element be- 
tween different (charge) ice configurations is then 
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(18) 



We will return to this model below in the context of predic- 
tions for experiment and simulations performed at finite tem- 
perature^. 

The manifold of ice configurations on the pyrochlore lattice 
is equivalent to the set of possible close-packed loop coverings 



of the diamond lattice^' . Exactly parallel arguments, leading 
to a formally identical Hamiltonian, can also be constructed 
for the closely related quantum dimer model on the diamond 
lattice^^i^. This model also exhibits a quantum U{1) liquid 
ground states for a smaller — but none the less finite — range 
of parameters 0.75(7 < M ^ 5 ll60ll6lll . 



B. Electromagnetism on a diamond lattice 



The mappings described in Section HTAl permit us to reduce 
complicated interactions between magnetic ions to a problem 
of tunnelling between spin configurations obeying the "ice 
rules" [cf. Fig. |2l. If we think of these spins as field lines 
of a fictitious magnetic field B, these rules can conveniently 
be written as 

V -B 

This naturally suggests an analogy with magnetostatics, with 
magnetic field lines constrained to lie on the bonds of a dia- 
mond lattice [cf. Fig.|6ta)]. And in the presence of tunnelling 
between ice configurations, this analogy can be extended to 
a fully dynamical quantum electromagnetism. Here we re- 
view the mapping from an ice with tunnelling, to a compact, 
U{1) lattice gauge theory, before moving on to an analysis 
of its "photon" excitations [Section HTCl and spin correlations 
[Section [ill CI . In so doing we follow closely the arguments 
of Hermele et atf,, but recast the discussion in terms of the 
magnetic fields B usually associated with the spins of a spin 
ice. 

We begin by transcribing the spin variables of "Htunneiiing 
[Eq. (fT2b l in terms of a quantum rotor variable 9i, and its con- 
jugate number operator rii 



1 

rij 

2 



where 



SJ = y^icxp [i9i\\/l - n. 
= Vl — n-i cxp [—i9i\\/n 



(19) 

(20) 
(21) 



(22) 



The number operator rii could equally be associated with the 
density of (hard-core) bosons in a charge ice, and in order to 
remain in the physical subspace where Ui — or 1, we add 
the term 



(23) 



to the Hamiltonian, subsequently taking the limit U ^ 00. 
With this restriction in place, the Hamiltonian becomes 



= I ^(n.- 1/2)2 

i 

O 



(24) 
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It is from this rotor form of the Hamihonian that we will make 
the passage to a gauge theory on the diamond lattice. 




curl < curl > curl = 



FIG. 6: (Color online). The different fields used in Section lllBl to 
construct lattice gauge theory of the spin liquid state, and the differ- 
ent lattices on which they are defined, (a) The "magnetic" field B^r' 
[Eq. (I26H, and its conjugate field CJ^r' [Eq- J27i l1 are defined on the 
links of the diamond lattice, shown here in red. Each link of this di- 
amond lattice corresponds to a site of the original pyrochlore lattice, 
and Byyi encodes the orientation of the spin on this site, (b) The com- 
pact (7(1) gauge field Ass' and the conjugate "electric" field £ggi are 
defined on the links of a second, dual, diamond lattice, shown here 
in blue. The midpoints of these bonds also form a second, dual, py- 
rochlore lattice, corresponding to the centres of hexagonal plaquettes 
in the original pyrochlore lattice, (c) An illustration of taking the 
lattice curl on the hexagonal plaquettes of the diamond lattice. The 
resulting vector lives on the links of the dual, diamond lattice. 

The site i of the pyrochlore lattice can be thought of as the 



midpoint of the bond r — ^ r' of a diamond lattice [cf. Fig.|3. 
Since this diamond lattice is bipartite, it is possible to define 
directed variables on these bonds 

Brr' = —Br'r Grr' ~ —Qr'r (25) 

through the mapping 

Brr' = ±(n,- (26) 

Grr' = ±0, (27) 

where the sign is taken to be positive if r belongs to the A- 
sublattice, and negative if r belongs to the S-sublattice. Tak- 
ing this convention into account, we are left with a pair of 
canonically conjugate variables 

[Qrr',Br"r"'] = i {Srr"5r'r"' ^ Srr"'Sr'r") ■ (28) 

The field, Brr' will take on the role of a magnetic field in 
our lattice field theory. However in order to recreate "electro- 
magnetism" we need also to discover an analogue to the elec- 
tric field. The missing field. Ess', inhabits the bonds s — s- s' of 
a second diamond lattice, interpenetrating the first [Fig.|6lb)]. 
It is defined through a lattice curl 

^^ss' = (Vo X g)ss' = J2 ^rr' (29) 

o 

where the sum is taken with anticlockwise sense around 
the hexagonal plaquette of pyrochlore lattice sites encircling 
the bond s ^ s'. It follows that £ss' is also a directed variable 

£ss' = -Ss's (30) 

We are now in a position to transcribe "Hrotor completely in 
terms of "electromagnetic" fields 

^rotor = ^ ^ - 2.9 ^ COS (f ,,0 (31) 

(rr') (ss') 

where the sum J2{rr') '^^^^ bonds of the original dia- 
mond lattice, while the sum X]{ss') '^^^^ bonds of the 
second, dual diamond lattice. The fact that the Hamiltonian is 
invariant under the transformation £ss' — > ^^ss' + 27r makes it 
evident that this theory is compact. It is also important to note 
that each of the components of the total magnetic field 

iB^,By,B,) = J2 Brr' err' (32) 
(rr'> 

where err' is a unit vector directed from r to r', is a conserved 
quantity under the dynamics of Hmtor [Eq. dSTTil. More gen- 
erally, reversing the sign of Brr' on a closed loop of spins will 
tunnel one ice configuration to another, without changing the 
total magnetisation of the system. 

In deriving Eq. (|3TI ). we have assumed that the ice rules 
hold, i.e. 

(V-6)r = X^'Srr' =0 (33) 
(r') 
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where the sum J2{r') '^^^^ ^^^^^ neighbouring r. This 
condition is automatically satisfied if we write Brr' as the lat- 
tice curl of a gauge field ^ss' ■ However we must also respect 
the requirement that the field B^r' take on half-integer values 
[cf. Eq. (|26])]. This can be accomplished by introducing a 
static background field B^^,, taken from any spin configura- 
tion which satisfies the ice rules, and writing 

{Brr'-B",,,) = (VoX^),^, (34) 

to give 

(rr') 

-25^cos(£ss') (35) 

(ss'> 

The fields £ss' and ^ss' are canonically conjugate 

[£ss', A"S"'] = i {Sss"Ss's"' - 5ss"'Ss's") (36) 

Moreover, the theory has a local gauge symmetry since one 
can make the transformation 

As' ^ ^ss' + As - As' (37) 

on any bond without changing the value of (Vo x A)rr' — 
each value of Ag occurs twice, with opposite signs. The situ- 
ation now bears more than a passing resemblance to quantum 
electromagnetism. 

At this point a subtlety enters the problem. In passing from 
Eq. (O to Eq. (l35l l. we have performed a series of changes of 
variable without making any new approximations. However 
it still remains to take the limit C/ — > oo. If the "magnetic" 
field Brr' were an integer variable, it could be eliminated from 
the problem by setting Brr' = on all bonds. This would 
be energetically favourable at large U/g, and would imply a 
phase transition from a spin liquid phase at small U/g , into a 
phase in which spinon excitations (magnetic monopoles) were 
confined at large U / g [cf . ItH] . However the fact that Brr' 
takes on half-integer values "frustrates" the lattice theory, and 
makes it possible for a spin liquid phase to survive in the limit 
U ^ oo. 

Keeping this in mind, we now follow Hermele et al^ in 
assuming that an average over fast fluctuations of the gauge 
field a) softens the restriction that Brr' take on half-integer 
values and, b) restricts fgs' to small values. Provided that both 
of these assumptions hold true, we can drop the reference field 
and expand the cosine in Eq. QST l, to obtain 

^u(i) = I E [(^o X A\A" + f E ^ss' (38) 

(rr') (ss') 

where both the normalisation of the field Brr', and the param- 
eters of 7iu{i) rnay be renormalized from their bare values 
\Brr'\ ^ l/2,U ^ U, JC ^ g. This, finally, is the Hamil- 
tonian for non-compact quantum electromagnetism on a dia- 
mond lattice. 



At first sight, the final step of this derivation might seem 
to involve an uncomfortably large leap of faithiS. However 
this will be justified a posteriori in Section IIIIBI by the ex- 
cellent, quantitative, agreement of the predictions of T-Lu(i) 
[Eq. ( l38b l with quantum Monte Carlo simulation of the mi- 
croscopic model T-L^ [Eq. (fT6]l1. In order to extend this com- 
parison to finite values of the control parameter /i, we will 
augment H (7(1) with a term 

-^^ud) = Y E [(^° ^ ^ -^))ss'] (39) 

(ss'> 

which mimics the effect of the "RK" potential [Eq. ©I- Since 
this term is permitted by the gauge symmetry, in principle it 
might also be generated dynamically by an average over fast 
fluctuations of Ass' ■ 



C. Constructing the photon 

The lattice gauge theory described in Section HTB] supports 
three types of excitation : magnetic charges (point sources 
of B) and electric charges (point sources of £), together with 
photons (transverse excitations of A) which mediate Coulomb 
interactions between these emergent charges ll46ll47ll . 

The magnetic charges are the magnetic monopoles of the 
classical theoryi^, now quantised and endowed with dynam- 
^^ ^40,74^75 ^ They correspond to the "spinon" excitations of the 
spin liquid. Since they involve spin configurations lying out- 
side the ice manifold, they have an energy gap 

[cf. Eq (fT0tl. The electric charges are gapped, topological 
excitations which can be constructed as a wave packet of ice 
configurations with suitably chosen phase a"^^'"^^ . These also 
have an energy gap 

Ae^JC^g^ 124/ 

[cf. Eq (O]. 

However the energy of the photons vanishes linearly at 
small wave vector 

a;(k^ 0) = c|k|. 

and being gapless, the photons will control the low energy 
and low temperature properties of the system. We therefore 
concentrate on exploring the consequences of the photons in 
this paper, leaving other excitations for future work. 

In what follows, we will explicitly construct a photon basis 
for the lattice gauge theory developed in Section IIIBI with a 
view to calculating the spin-spin correlation functions of the 
original model of a quantum spin ice. We take as a starting 
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point 



^U(l) = Y E [(^O -^)(r,n 



1 2 



2/C 



sG A' ,m 



+ Y E [(^ox VoX^)(,,„J (40) 

where we have used the fact that, in the absence of electric 
charges 

1 dA(^s,m) 



(s,m) 



(41) 



/C dt 

To avoid double counting of bonds, the sums over diamond 
lattice sites {r} and {s} are restricted to a single sublattice, 
with bonds labelled 



(r,n) = (r,r + e„) , (s, to) (s, s + e,„) 



where 



eo 
ei 

62 



(1,1,1) 



(-1,1,-1) 



ao 



63 = ^(-1,-1,1) 



(42) 



and ao is the linear dimension of the cubic unit cell of the 
lattice, shown in Fig.|5] 

We proceed to quantise ^s,„ by analogy with conventional 
electromagnetism, introducing a Bose operator 

[a^,a{,] = Sxx' 

where the four sites of the tetrahedron in the primitive unit cell 
of the pyrochlore lattice translate into four bands A = 1 ... 4. 
We write 



A 



(s,m) 



EE 



^A(k) 



k A=l 

X (exp [-ik ■ (s + e„j/2)] 77„A(k)a^(k) 

+ exp [ik ■ (s + e^/2)] rjl^{k)a{{k) 

(43) 

where the sum X]a=i ^^^'^ ^'^^^ branches of photons 
and ?2(k) is a unitary, 4x4 matrix whose columns, 'i]\{k), 

play the same role as the polarisation vector in conventional 
electromagnetism. By obvious extension 
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(s.m) 



EE 



k A=l 



c^A(k) 



(exp [-ik ■ (s + e„i/2)] 77„iA(k)a^(k) 
- cxp [ik • (s + ern/2)] r/*„ (k)at (k) 



(44) 



The Hamiltionian (Eq.l40b is already quadratic in a^. What 
remains is to eliminate all terms which do not conserve pho- 
ton number, by constructing a suitable matrix /^^^(k). To do 
this, we need to evaluate the Fourier transform of the lattice 
curl (Vo X „■). This operator is defined on a six-bond 
plaquette, composed of pairs of bonds which enter with op- 
posite signs in the directed sum around the plaquette [Fig. |2] 
Fig.|6l. These bonds have midpoints located at 



r - e„/2± h„ 



where 



Hence 



ao e„ X e„ 



V8 



(45) 



(Vo X 



/V E E 



IC 



k A=l V ^ ^ 



X < cxp[— ik • (r — e„/2)]a^(k) 

X ^(-2isin(k ))?/mA(k) 

ra 

+ exp [ik ■ (r — e„/2)] a^(k) 
x^(2isin(k.h„™))?7^„(k)| (46) 

where, by inspection, h„„ = 0. 

We can rewrite the sum in Eq. ( |46b in a more conve- 
nient form by introducing an Hermitian, anti-symmetric ma- 
trix 



'2i X 



/ sin(k • hoi) sin(k • ho2) siii(k • hos) 

— sin(k-hoi) sin(k • hi2) siii(k • his) 

-sin(k-ho2) -sin(k-hi2) sin(k • h23) 

V-sin(k-ho3) -sin(k-hi3) -sin(k-h23) 



(47) 



acting on the four component vectors 77a (k). 

Since Z.ik) is Hermitian, we are free to construct the matrix 

2(k) from the eigenvectors of Zjk), such that 



^(k) 




= CA(k) 




(48) 



A specific choice of ^(k) corresponds to a choice of gauge, 

since using Eq. ( 1431 ). the divergence of ^ss' is now fixed. The 
choice here, which is made for maximum convenience in con- 
structing the photon dispersion, is the radiation (or Coulomb) 
gauge 



(49) 
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It follows from Eqs. ( |46] ) and ( |48] ) that 



k A=l V ^ 

X I exp[-ik • (r - e„/2)]aA(k)CA(k)7?„A(k) 

+ exp[jk • (r - e„/2)]al(k)CA(k)77L(k)^ . (50) 
Squaring and summing over r and n, we arrive at 



(r,n) 



k A=1A' = 1 



XW-^J— ^CA(k)CA'(k') 



CLiA(k) Y a;A'(k) 
X |aA(k)aA' (-k)CA (k)CA' (-k) ^ Vnx{k)vnX' (-k) 
+a{{k)al (k)CA(k)CA' (-k) J2 riU^Wx'ni-^) 

■>l 

+aA(k)al, (k)CA(k)CA' (k) ^ r;„A(k)r;I,„(k) 

n 

+al(k)aA' (k)CA(k)CA' (k) J2 Vln{^)VnX' (k) j . (51) 



This rather dense expression can be simpUfied using the uni- 
tarity of ^(k) 



Y,VU^)VnX'{k)^6. 



_^ .... AA'- 
n 

and the fact that 

K-k) = ^(k)* 
from which it follows that 

Uk) = CA(-k). 

Whence, 

E(Vox <„., = ! EE 

(r,n) k A=l ^ 

x(aA(k)aA(-k)+al(k)al(-k) 



(52) 



(53) 
(54) 



+4(k)aA(k)+aA(k)al(k) 



(55) 



Applying the same procedure again to Eq. (l50l ). we find 
E(VoxV„x^)U,4EE^ 



(s,m) 

x(aA(k)aA(-k)+at(k)4(-k) 



+4(k)aA(k)+aA(k)at(k) 



(56) 



The remaining, electric field, term in Hq [Eq. (|38] )1 yields 



fdA 



(s,m) 



-^EE-^(k) 



k A=l 



x I - aA(k)aA(-k) - al(k)4(-k) 

+4(k)aA(k)+aA(k)al(k)|. (57) 
Inserting all of this into the Hamiltonian Eq. ( l40l i gives 



'^u(i) "EE 



fWCCxiX)^ >V/CCA(k)^ 
^ 4L.A(k) ^ 4L.A(k) ^ 4 



k A=l 

x (aA(k)a^(k) +at(k)aA(k 

/ Z^^CA(k)^ >V/CCA(k)4 _ 
\ 4c^A(k) 4c^A(k) 4 

aA(k)aA(-k)+at(k)at(_k) 

To diagonalize the Hamiltonian we require 

Z^/CCA(k)2 _ W/CCA(k)4 c^A(k) 



4wA(k) 
which implies 



4wA(k) 



4 / -, 

^u(i) = E E ^^(1^) 4(k)aA(k) + - 



k A=l 

with dispersion relation fixed by Eq. 



C.A(k)=/C^/|CA(k)2 + ^CA(k)4. 



(58) 



(59) 



(60) 



(61) 



All that now remains is to determine the eigenvalues of the 
matrix Z{k), CxO^)- We find 



Ci(k) = +V2 J2sm{k-h^^y 

y mn 

C2(k) = -\/2 /^sin (k-h™„)' 



C3(k) 

C4(k) 



(62) 

(63) 

(64) 
(65) 



It follows that the four bands of excitations (\ (k) correspond 
to two, degenerate, physical photon modes, and two unphysi- 
cal, zero energy modes. The unphysical modes arise because 
of the gauge redundancy in A and make no contribution to 
either the Hamiltonian or to any gauge invariant correlation 
functions. 

Keeping only the physical photon modes from Eq. ( |60] |, we 
finally arrive at 



^U(i)=EE^(k) (a{{k)a,{k) + l 



k A=l 



(66) 
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FIG. 7: (Color online). Dispersion a;(k) of magnetic photon exci- 
tations, calculated for the lattice field theory Eq. ( 140) in the "quan- 
tum ice" limit W/K. — >■ 0. The dispersion is plotted in the (ft, h, I) 
plane, following Eq. l l67t . The dispersion is linear in |k| in the long- 
wavelength limit, with a speed of light c = VWCao- 



where A now has the interpretation of the polarisation of the 
photon. The photon dispersion c<;(k) is independent of polari- 
sation and can be written 



a.(k) = /C\/^C(k) + ^C(k)2 



(67) 



where 



C(k) = Cl(k) - -C2(k) = V2 /^Sin(k • hrnnf (68) 

y nin 

with h,„„ defined by Eq. (05]). 

For all U/K, > the photon dispersion is linear in the 
long-wavelength limit 

LuCkM 0) M VWCao\k\ (69) 

This means that there is a well-defined speed of light 

c=VWCaoh~^ (70) 

where we have restored the dimensional factor of h. 

However in the limiting case U/K, — 0, c — > 0, and 
the dispersion of the photon becomes quadratic in the long- 
wavelength limit 



cj(k) « VWlCal\k\- 



(71) 



Precisely this limit is realised at the RK point /^i = g of the 
quantum ice model [Eq. (fTSIll. and defines the boundary 
of the quantum C/(l) liquid phase^^i^. The photon dispersion 
relations in the two extreme cases U /K, = and W//C = 
are plotted in Fig.|2]and Fig. [8] 



D. From photons to structure factors 

Spin correlations in real materials can be measured directly 
by neutron scattering. Here we convert the analysis of pho- 
tons in Section HTCl into concrete predictions for the dynami- 
cal structure factors measured in such an experiment. We also 



consider the structure factors which might be measured in, 
e.g.. X-ray scattering experiments on a charge ice of the type 
considered by Banerjee et al.—. Specifically, we will consider 

= Idie— *(S"(-k,OS^(k,0)) (72) 



and 



5charge(k,w) j dte-'^\n{^\i,t)n{\i,Q)) (73) 

Possessing the full photon wave function [Eq. ( l43] ll permits 
us to calculate these dynamical structure factors on a lattice, 
passing directly from the correlations of A{^s,m) to those of 

S(r,n) Orn(r_„). 

We first consider the charge ice and, following |59l intro- 
duce an additional (dimensionless) scale factor k < 1 to take 
account of any renormalization of the field B when an aver- 
age is taken over fast fluctuations of A^s m) [cf- Eq. ( [35] l to 
Eq. dUl)] 

S'charge(k,t) = ^ (6„ (-k, t)^,„ (k, 0)) (74) 



where 



^^^^^^Y. cxp[-«k • (r + e„/2)]S„(r, t) (75) 



with 



S„(r) = B(r-e„/2) = (VoX^)(r_„) 



The time evolution of A^s.m) follows directly from "Hyj-j^-j 
[Eq. (|66l)] 



:CA(k) 



4 V 
(77„A(-k)aA(-k)e-'"-('')* 

+77L(k)al(k)e-^('^)*) . (76) 



such that 



charge 



mn A 



rt'^A(k) 



x(a;,(k)4(k)+al(k)a^(k)) 
y.5{u-ujx{\<i)) {11) 

where we have dropped all terms which fail to preserve photon 
number or polarisation. 

The unphysical photon polarisations A = 3,4 do not con- 
tribute to Eq. ( TTTI i. since 



CA(k)V^A(k)|A=3,4=0 

For the physical polarisations A = 1,2 



(4(k)a^(k)) = 



1 



e T — 1 



= nB(k) 



(78) 



(79) 
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FIG. 8: (Color online). Dispersion cj(k) of magnetic photon ex- 
citations, calculated for the lattice field theory Eq. ( 140) in the limit 
U/K. — >■ 0. The dispersion is plotted in the (h, h, I) plane, following 
Eq. ( 167b . The dispersion is quadratic in k the long wavelength limit. 
This situation is realised in the microscopic "quantum ice" model 
[Eq. ([T6ll], at the RK point, fi ^ g. 



since the photons are bosons. Noting that 



It is this form of the result, evaluated at the discrete set of wave 
vectors {k} appropriate for a finite-size cluster with given 
boundary conditions, which we will fit to simulation results 
in Section lmB] 

Calculating the dynamical structure factor 5^f^(k, lu) for a 
spin ice means generalising Eq. (ISTI i to take account of neu- 
tron polarisation, and the local easy axes of spins in a spin 
ica^. However the underlying field-theoretical description of 
the problem Hy^-^j [Eq. ( l40b l is unchanged, and the two re- 
sults differ only in the way in which the contraction of fields 
on different sublattices (S„(— k)S„j(k)) contribute to corre- 
lation functions. In a charge ice we simply sum over m, n as 
in Eq. ( 1741 ). In a spin ice we must account for the easy axes 
which lie along the vectors e„ [Eq. ( l42l and then calculate the 
projection of the spin along the axis of interest^. Thus, the 
equal time structure factor is 



e™-a e„-/3 (e„(-k)B„(k)) (86) 



^CA(k)2ry„a(k)7yL=^(^(k)2), 



(80) 



we arrive at a result for the dynamical structure factor of a 
quantum charge ice 



•Schargelk, Uj) 



2 a;(k) 



^ ^ sin (k • h„;) sin (k • h„;) 



mn I 

+ Siuj + ujx(k.))nB(k) 



(81) 



where the vectors h„m are defined by Eq. (1431 ). 

In comparing with quantum Monte Carlo simulation, we 
will also make extensive use of the zero-temperature, equal- 
time (i.e. energy-integrated) structure factor 

•Schargelk, t = 0)t=0 ^ j -Schargelk, w)t=0 (82) 

This can be written as a function of just two, dimensionless, 
ratios of parameters, U and W 



'S'chargeCk, i = 0)t=0 



^o(k) 



(83) 



yZ^C(k)2 + WC(k)4 
where C(k) is defined by Eq. ( |68] |. 

^0 (k) = ^ ^ sin (k • h™/ ) sin (k • h„z ) , (84) 



mn I 



and the dimensionless ratios of parameters are given by 
— U — W 

U^—r , W 



ICk* 



(85) 



where a and /? are unit vectors in the a and j3 directions. 

Following the same procedure as described above for the 
charge ice we come to the general result for the dynamical 
structure factor in a spin ice 



2 uj(k) 



sin (k • h^;) sin (k • h„;) 



mn / 



X • a e„ • /3 U - ujxik)) (1 + nB(k)) 



+ 6{uj + uJxCk))nB{'k) 



(87) 



For concreteness, where we come to plot results, we will 
follow the conventions of Fennell et al.—, who used neutrons 
with polarisation parallel to 

= (1,-1,0) 

to measure the energy-integrated structure factor 
S"p^^(k,t = 0), for transfered momentum k in the {h,h,l) 
plane. We also follow the conventions of Ref. |^ in choosing a 
coordinate system in which 



X II k , y II n,. X k , z II rii., 



(88) 



and consider the "spin-flip" channel S'^p^|^(k, oj). In this con- 
vention, the non spin-flip channel measures S^p^^^Ck, u). 

It follows from Eq. dSTb that the dynamical structure factor 
in the spin-flip channel is given by 



^2 ^ 



2 w(k) 



y ^ y ^ sin (k • h,„;) sin (k • h„;) 



mn I 

Bm • (n,, X k) \ fen - (n^ x k) 



|(n, xk)| ; V IKxk)| 
x{S{uj- ujxm (1 + nsik)) +Siuj + ujx{k)) nsik)) 
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and the corresponding zero-temperature, energy-integrated 
(i.e. equal-time) structure factor is 



2 Lj{k) 

X ^ ^ sin (k • h„,7) sin (k • h„i) 

nin I 

e„j • (n^ X k) \ /e„ • (n^ x k) 



X 



fn„ X k) 



(n,^ X k) 



(89) 



Once again, we will make extensive use of this result when 
comparing with quantum Monte Carlo simulation. At finite 
temperature we obtain, for the energy integrated structure fac- 
tor in SF channel 

3"^ (kt-0)- — ^ coth ( ^ 
X ^ ^ sin (k • h,„/) sin (k • h„/) 

^' e„ ■ (n^ X k) ^| |^ e„ ■ (n^ x k) '^ ^^^^ 



(n,y X k) 



(n^ X k)| 



Where neutron scattering is performed with unpolarized 
neutrons, experiments measure an average over different com- 
ponents of the dynamical structure factor 



/(k,a;) PC ( 



'-j^)s:iiK.) (91) 



This is the result plotted where we illustrate photon disper- 
sions in Fig. [3] and Fig. [T3] The corresponding quasi-elastic 
(energy integrated) form of Eq. (|9T1 ) is given by 



ct/3 



)5rpC(k,t = 0). (92) 



It is important to note the predictions for quasi- 
elastic neutron scattering S^^-^^ik.ijj) [Eq. ( [89l )1 and 
/(k, w) [Eq. (|9T1)1. and their energy-integrated counter- 
parts Sl^-^lk^t = Q) ^] and /(k) [Eq. ^] only 
include contributions from the low-energy photon excita- 
tions of a quantum spin ice. In a real quantum spin ice 
material, their would also be contributions at higher energy 
from gapped "electric charges", and magnetic monopole 
excitations. These are not treated in the present theory. 



quantum ice model Ti,^ [Eq. (fT6]l1, by making explicit com- 
parison with the results of zero-temperature quantum Monte 
Carlo simulation. 

However before considering results on a lattice, it is useful 
to ask how correlations in quantum spin ice might differ from 
those in a classical spin ice, within a simple continuum field 
theory. This is considered in Section IIII Al We then turn to 
simulation of the lattice model T-L^ [Eq. ( fTSI ll in Section lTllBI 
demonstrating that the lattice field theory 'Hu(i) [Eq. ( I38] l1 
provides an excellent quantitative description of the results 
for Cpfnlk,^ = 0). In Section |IIIC] we use the same lattice 
field theory to make predictions for the magnetic photon ex- 
citations which could be observed in inelastic neutron scatter- 
ing experiments. Finally in Section lTlI Dl we use the finite-size 
scaling of ground state energies in simulation to put an abso- 
lute scale on the speed of light c associated with these mag- 
netic photons. Throughout this analysis we set h = \, restor- 
ing dimensional factors of h only where we quote results for 
the speed of light. 



A. Structure factors within continuum theory 

The long-wavelength properties of a quantum C/(l) liquid 
are well-described by a continuum field theory of the form 
considered in Ref . |4g 



5eff = -!- / dtcPr 

OTT 



£{vf - c^Bivf 



Pc V X B{v) 



(93) 



This therefore provides a convenient starting point for dis- 
cussing the evolution of spin correlations in quantum spin ice. 
We emphasise that such a theory can be derived as a con- 
tinuum limit of Hy^j [Eq. (l40bl^^. And where we go on to 
make comparison with quantum Monte Carlo simulation in 
Section lTllBI we will use the appropriate results on a lattice, 
i.e. Eq. ^ and Eq. 

For pc = 0, iSeff reduces to the familiar Maxwell action 
of quantum electromagnetism. Crucially, this action supports 
photon excitations with dispersion ci;(k) = c|k|. The addi- 
tional term (V x B{v))^ is invariant under gauge transfor- 
mations .4(r) — >■ A{t) + V(/i(r), and is an irrelevant pertur- 
bation in the RG sensaii. However it introduces a new length 
scale into the problem 



(94) 



III. "ELECTROMAGNETISM" IN A QUANTUM SPIN ICE 

ATT = 



The arguments presented in Section HTBl explain how a spin 
liquid state with correlations described by an effective electro- 
magnetism can arise in a quantum spin ice, but stop short of 
offering proof that this happens in any real material or micro- 
scopic model. In what follows, we validate our use of Gaus- 
sian electromagnetism Hu(i) [Eq. ( l38T l1 as a description of the 



which controls the curvature of the photon dispersion 




(95) 



and has an important impact on how correlations evolve as a 
function of distance. 

The role of Ac can most easily be understood in the limit 
c 0, where correlations of B{y) are controlled entirely 
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FIG. 9: (Color online). Comparison between the predictions of the lattice field theory [Eq. | |40H and quantum Monte Carlo simulation 

of the microscopic model "H^ [Eq. l|16H, for a quantum charge ice at T = 0. First column : equal-time structure factor Scharge(k, i = 0) 
calculated using Green's function Monte Carlo (GFMC) simulation of a 2000-site cubic cluster, for a range of /i ranging from n/g = 1 (RK 
point) to fi/g = (quantum ice). Second column : best fit of the finite-size (FS) prediction of the lattice field theory to simulation, following 
Eq. ( I83t . There is excellent, quantitative, agreement between theory and simulation for all values of fi/g. Third column : prediction of lattice 
field theory in the thermodynamic limit, for parameters obtained from fits to simulation. 
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by pc- Precisely this limit is realised in the microscopic model 
[Eq. ( fT6b l at the exactly soluble "RK" point fi — g. At the 
RK point, all ice configurations are degenerate, and the pho- 
tons have dispersion cj(k) = y/p^|kp ll46i47ll . Correlations 
of the magnetic field 

Cf^Jk) = (6,,(-k)6,(k)) (96) 

can be calculated from Eq. ( |93] ). and for c = these behave 
as 



8^ 



(97) 



exhibiting the pinch-point singularities characteristic of the 
"Coulombic", classical U{1) liquid phaseiMl. On Fourier 
transform, Eq. (|97] ) corresponds to dipolar correlations in a 
three-dimensional space 



(98) 



The quantum U{1) liquid phase, with its Unearly dispersing 
photons, is stabilised by the emergence of finite value of the 
speed of light c for fJ. < g |46,47 58,]. In this case, we find 



1 



(99) 



[cf. I47II48I1 . For wavelengths A < Ac, Eq. (|99]l reduces to 
Eq. ( |97l ), and the system exhibits "classical" dipolar correla- 
tions of the form Eq. (|98l ). However for long wavelengths 
A ^ Ac the additional factor of k in the numerator of Eq. (|99] | 
"hollows out" the pinch point singularities. In this limit, 
r ^ Ac, Eq. ( |99l ) corresponds to dipolar correlations in ?l four- 
dimensional space 



(100) 



the additional dimension arising because of fluctuations in 



timi 



,47.48 



We therefore associate Ac with the length-scale over 
which the system crosses over from "classical" ice correla- 
tions, decaying as 1/?''^, to "quantum" ice correlations decay- 
ing as l/r"*. 

The length-scale Ac will also play an important role where 
we compare the predictions of field theory with simulation of 
microscopic model T-L^ [Eq. (fT6] l1 as a function of We can 
gain some insight into the /i dependence of Ac, from degener- 
ate perturbation theory about the RK point^i^^iffi^. We find 
that ~ {g — /i), while pc ~ const., and it follows from 
Eq. ( |94l i that Ac diverges as Ac ^ 1/ ^Jg — fJ,- Exactly at the 
RK point, where g = /i, Ac is infinite and correlations have 
the classical form Eq. ( l98T l at all length scales, as expected. 
However, as we move away from the RK point into the quan- 
tum liquid phase for p/g < 1, there will be a progressive 
evolution of correlations from classical (pinch points) at short 
distances to quantum (no pinch points) at long distances. This 
expectation is born out by quantum Monte Carlo simulations, 
described below. 



K 0.5 




FIG. 10: Parameters of lattice field theory "H^i) [Eq. as a 
function of n/g, from comparison with quantum Monte Carlo simu- 
lation of spin correlations in the microscopic model [Eq. ( I16b 1 
at T — 0. The dimensionless combinations of parameters 
U and W [Eq. J85H were obtained by fitting the predictions of the 
lattice field theory S'charge(k,f = 0)t=o [Eq. {SS}] to simulation re- 
sults at fixed n/g [cf. Fig[9). 



For the purposes of these simulations, Ac also sets the min- 
imum size of cluster which is needed to capture quantum ef- 
fects at a given /i. At /i = we find that Ac ~ 0.8ao, and 
hence a cluster of linear dimension L = 5ao {N = 2000) is 
comfortably big enough to observe the quantum spin Uquid 
phased. 



B. Comparison with Quantum Monte Carlo simulation 

We now turn to zero-temperature quantum Monte Carlo 
simulation of the microscopic model [Eq. (fT6lll. We have 
previously argued that this model supports a quantum U{1) 
liquid ground state for a range of parameters —O.bg < /i < g 
— cf. Fig. Hlandfsil In this earlier work, evidence for the 
ground state phase diagram was taken from the finite-size 
scaling of energy spectra. Our main tool here will be the equal 
time structure factor S'(k, t ~ 0), calculated from simulation, 
and from the lattice field theory Hyj^-^j [Eq. (|40b]. These two 
independent calculations are found to be in excellent, quanti- 
tative agreement, confirming the conclusions of IsH Making 
a direct comparison between the field theory and simulation 
also serves to put the field theory on a quantitative footing, 
providing information about the evolution of the parameters 
of the field theory as a function of the microscopic parameter 
/i. 

Simulations were performed using a Green's Function 
Monte Carlo (GFMC) technique based on the statistical sam- 
pling of ice configurations. This sampling is weighted using a 
variational estimate of the ground state wave function, which 
is optimised in a separate variational Monte Carlo (VMC) cal- 
culation. In this sense, GMFC can be thought of a systematic 
method of improving upon variational calculations. There 
is no sign problem associated with "H^, since all of its off- 
diagonal matrix elements are equal to or —g, with g > 0. 
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FIG. 11: (Color online). Comparison between the predictions of the lattice field theory 'Hu(i) [Eq. ( I40H and quantum Monte Carlo simulation 
of the microscopic model "Hp [Eq. ||16H, for a quantum spin ice at T = 0. First column : equal-time structure factor S^^^^{k,t — 0), as 
measured in neutron scattering by Fennell et al^, calculated using Green's function Monte Carlo (GFMC) simulation of a 2000-site cubic 
cluster for parameters ranging from /i/g = 1 (RK point) to fi/g — (quantum ice). Second column : best fit of the finite-size (FS) prediction 
of the lattice field theory to simulation, following Eq. i90\ . There is excellent, quantitative, agreement between theory and simulation for all 
values of fi/g. Third column : prediction of lattice field theory in the thermodynamic limit, for parameters obtained from fits to simulation. 
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Where simulations converge, the results obtained are numer- 
ically exact. Our implementation of VMC and GFMC cal- 
culations for quantum ice^^ exactly parallels our earlier work 
on the quantum dimer model on a diamond lattic o^°'^' , with 
correlation functions calculated using techniques described in 
Ref. We refer the interested reader to these papers for 
further details of the method. 

In the left-hand column of Fig. |9l we present GFMC sim- 
ulation results for the equal-time correlations in a quantum 
charge ice 

S'charge(k,t = 0)t=o = (n(k)n(-k))T=o 

Simulations were performed for a 2000-site cubic cluster 
possessing the full symmetry of the lattice, for parameters 
n/g = 1, 0.75, 0.5, 0.25, 0. 

The classical, dipolar correlations at the RK point 
\xl g = 1 are clearly visible as sharp "bow-tie" motifs in 
'S'charge(k, t ~ 0), Centred on pinch-points at k = (1, 1, 1), etc. 
As expected, these pinch points are progressively eliminated 
as — > 0, and quantum effects come to dominate the long 
lengthscale physics of the problem. This erosion of the pinch 
points is accompanied by a gradual redistribution of spectral 
weight, with high intensity regions evolving from a triangular 
into an oval shape. 

In the central column of Fig. |9l we present the best fit to 
simulation results obtained from the lattice field theory. Fits 
were made using the result S'charge(k, t = 0)t=o [Eq. (|83])1. 
evaluated for the same 2000-site cluster, as a function of the 
two dimensionless parameters U and W [Eq. (I85]l1. The two 
results are indistinguishable by eye, and differ maximally by 
a few percent, for values of k close to the Brillouin zone 
boundary. The quality of these fits implies that they can be 
used to accurately parameterize the lattice field theory "Hy^j 

[Eq. (|40]|1. and the values of U and W obtained are shown in 
Fig. [lO] We note that the values obtained at the RK point, 
U = Q and W = 1, are uniquely determined by the known 
form of correlations within the classical ice states^. A sepa- 
rate evaluation of the speed of light c oc \JlAK from finite size 
scaling of the ground state energy is given in Section lTlIDI be- 
low. 

In Fig. [TT| we show equivalent results for the equal-time 
structure factor of a spin ice 

5j;^,(k,i = 0)t=o = (Snk)Sn-k))T=o 

in the spin-flip channel considered by Fennell et alX. Super- 
ficially, these results look very different to those presented in 
Fig.|9] This is because the local easy axis is different for each 
of the four sublattices, leading to a staggering of correlations 
not present in the charge ice problem. However the informa- 
tion content of the two structure factors is exactly the same. 

At the RK point [ij g ~ 1, correlations are classical, and 
S'J'p^^(k, t = 0) exhibits a characteristic "snow flake" motif in 
the (/i, /i, I) plane, also seen in neutron scattering experiments 
on Ho2Ti207 i^. Pinch point singularities are clearly visible 
at the reciprocal lattice vectors k = (1,1,1), etc. 

Once again, these pinch points are progressively eroded as 
the system is tuned away from the RK point into the quantum 



tu(k)/(c a„ ') 




FIG. 12: (Color online). Relationship between the dispersion of 
the magnetic photon excitation cj(k) [Eq. | |67H . and the equal time 
structure factor S'fp^(k,t = 0) [Eq. ^\ in a quantum spin ice. 
The photon dispersion a;(k) in the (/i, h, V) plane is plotted above 
the corresponding equal-time structure factor, demonstrating how the 
photon disperses out of the (suppressed) pinch points at reciprocal 
lattice vectors. Note that the intensity of the scattering S'Jpy„(k, t = 
0) where a;(k) [Eq. i fTOST ll. Results were calculated within 
the lattice field theory [Eq. ( I40H for W = 0, with energy measured 
in units such that h= \. 

spin-liquid regime for [ij g < 1. Probably the most striking 
change, however, occurs at k = (0, 0, 0). Here, for a classical 
spin ice 

SspJ^ -^0,t = 0)t=o ^ const. 
However, in a quantum spin ice, 

and spectral weight is progressively excavated from the region 
of reciprocal space around k — (0,0,0) for n/g < 1. This 
has important consequences for the evolution of correlations 
at finite temperature, discussed in Section IIV Al and for the 
uniform magnetic susceptibility, discussed in Section ITV CI 

We wish to emphasise that the results shown in Fig.fTTIare 
not the outcome of separate simulations of a quantum spin ice. 
They are taken from the same simulations of the quantum ice 
model Hfj, [Eq. (fTSbl, recast in the coordinates appropriate for 
a spin ice. It follows that the parameters obtained from fits to 
field theory at finite size are exactly the same as those for a 
charge ice, given in Fig.fTOl 

C. Seeing the light : photons and inelastic neutron scattering 

Inelastic neutron scattering provides a direct method of 
measuring the dynamical structure factor S'^f^(k,aj), and so 
of resolving photon excitations in a quantum spin ice. These 
photons disperse linearly out of those reciprocal lattice vec- 
tors where pinch points are observed in quasi-elastic scatter- 
ing experiments. However, since these experiments measure 
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i = (-0.5,-0.5 -0.5) r = (0,0,0) r = (2,2,2) r = (2,2.0) (1,1,0) r = (l,l,l) 

k 

FIG. 13: (Color online). Ghostly magnetic "photon" excitation as it might appear in an inelastic neutron scattering experiment on a quantum 
spin ice realising a quantum ice ground state, for a series of cuts along high symmetry directions in reciprocal space. The prediction of the 
lattice field theory [Eq. J40H for inelastic scattering by unpolarized neutrons, /(k, uu) [Eq. l |91| l1 has been convoluted with a Gaussian of 

variance 0.3 c Oq ^ to represent the finite energy resolution of the instrument. The intensity of scattering vanishes as — !> 0, and is strongest 
at high energies. Energy is measured in units such that ft = 1, and the photon dispersion calculated for W = 0. 




(h, h,0) 

FIG. 14: (Color online). Prediction of the lattice field theory 
'Hy(i) [Eq. for quasi-elastic neutron scattering performed us- 
ing unpolarised neutrons, for comparison with Fig. [13] Results for 
/(k) are taken from Eq. ( 192) . and calculated for W = 0. The path 
within the [h, h, I] plane used for plotting the photon dispersion in 
Fig.[T3]is shown using unbroken black arrows, with Brillouin zone 
boundaries marked as dashed white lines. 



the energy integral of the dynamical structure factor, the sup- 
pression of pinch points in a quantum spin ice at T = has 
important implications for the observation of its photon ex- 
citations. Specifically, for non-interacting photons, the sup- 
pression of energy-integrated structure factor must imply the 



suppression of the weight in the photon peak itself. This is 
illustrated in Fig. [12] 

To see how this works, we consider the result for 
the dynamical structure factor in a quantum spin ice 
5^f^(k, a;) [Eq. (|87])1. in the (physically relevant) limit where 
W = 0. In this case weight in the photon peak is determined 
by the ratio 

c.(k) 

where, 

S'o'^(k) = ^^sin(k • h„;)sin(k • h„;) 

mn / 

X (^e,„ • (^e„ • (101) 

and 

ujxik) = VWCCxik). (102) 

We can use the spectral representation of Z.(k) [Eq. (IHOt l to 
write 

sin(k • h„;)sin (k • h„;) 

Since the only contributions to the RHS of Eq. ( 11031 ) come 
from the two dispersing modes A = 1, 2, [cf. Eq. (|78])1. 
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Eq. ( llOlb simplifies to 



4 /CZY 

A— 1 mn 

X I e™ • a ) I e„ • /3 



(104) 



Expanding in the first Brillouin zone, for k w 0, we find 

E'7mA(k)77j„(k) (e,n ■ aj ( e„ • ) « ^ 



for a = f3 = y, z and zero otherwise. It follows that 
oc w(k) (5(a; — i^(k)). 



(105) 



Therefore at low energies, in the first Brillouin zone, inelas- 
tic neutron scattering experiments will resolve the magnetic 
photon excitation as a ghostly, linearly dispersing peak, with 
intensity vanishing as / cx ixi(k), as noted in IITOll . However 
at higher energies and in other Brillouin zones, the momen- 
tum dependence of 77mA(k)?7^^(k) in Eq. ( 1104b will lead to a 
significant variation in the intensity of the signal at fixed lu. 
This behaviour is illustrated in Fig. [13] where we have plotted 
the intensity of scattering /(k, w) [Eq. ( I9l1 i1 for an experi- 
ment performed using unpolarised neutrons. The correspond- 
ing quasi elastic scattering, and the path within the [h, h, I] 
plane, are shown in Fig. ( fT4l i. 

The phenomenology of this photon excitation stands 
in stark contrast to conventional antiferromagnets, whose 
linearly-dispersing spin-wave excitations have the greatest in- 
tensity approaching the zero-energy magnetic Bragg peak 
associated with magnetic order The difference between 
these two problems stems from the fact that the photon is 
a quantised excitation of A, while neutron scattering mea- 
sures correlations of B. The lattice curl needed to relate 
one to the other introduces additional factors of C\0^) in 
S'"pf„(k, uj) [Eq. (HZl)], which leads to the suppression of spec- 
tral weight at low energies. 

A better point of comparison, in fact, is the scattering of 
neutrons by real photons in models of a hot early universe^. 
In both cases, photons are associated with a periodically fluc- 
tuating magnetic field, transverse to the direction of their 
propagation. And in both cases neutrons scatter inelastically 
from these locally fluctuating magnetic fields. In a spin ice, 
this scattering can occur in both the spin-flip (SF) channel, 
in which case there is a transfer of angular momentum to the 
sample, and in the non spin-flip (NSF) channel (cf. Fig. [131 
Fig.[T4land Fig.[T6]). It is also interesting to note that the same 
phenomenology of linearly-dispersing excitations, with a van- 
ishing spectral weight at long wavelengths, is encountered in 
quantum spin nematics^**. In this case, low-energy spin fluc- 
tuations are controlled by a time derivative of the underlying 
nematic order parameter—, and so vanish for ut 0. 



D. Estimating the speed of light 

The signal feature of the quantum U{1) liquid is its photon 
excitations. One important consequence of these, so far as the 
simulation of finite-size systems is concerned, is a character- 
istic finite-size correction to the ground state energy per site 
Eq/N, coming from the zero-point energy of the photons 



N 



1 

N 



[Eo{L) - Eoi<x)] ^ x,L-^ + 



(106) 



where L ~ 3 is the linear dimension of the cluster, and 
the coefficient xi is proportional to the speed of light c [cf. 
I61I1 . This means that it is possible to extract the speed of light 
from the finite-size scaling of the ground state energy found 
in simulations of H^j^ [Eq. ([T6bl, shown in Fig. [15] 

Approaching this problem from the lattice field theory 
'H'u(i) [Eq. ([40])], we know that 



(107) 



where the dimensionless parameter U = can be deter- 
mined separately from fits to structure factors (cf Fig. [Tot . We 
also have enough information from the fits to the structure fac- 
tor to evaluate the sum 
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(108) 



where is the ground state energy per site found from Monte 
Carlo simulations. 

For = the LHS of Eq. dlOSI ) does not depend on L. 
This is consistent with simulations of the microscopic model 
9- 

We write 



atp~g. For U > Owe expect a scaUng law ~ for large L. 



and it follows that 



with 
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(109) 



(110) 



(111) 



where the coefficients xi and X2 can be found from the finite- 
size scaling of the ground state energy in simulation [Fig. [T5]| , 
and through the numerical evaluation of the in Eq. (11081 ) 
for a finite-size system. 

We find that, for < /i < 1, the evolution of the speed of 
light as a function of /i is well-described by 



where 



Sfi = l-fj./g 
a = 0.22 5^ al 
p = 0.13 g^al 



(112) 

(113) 
(114) 
(115) 
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FIG. 15: (Color online). Finite-size scaling of the finite-size cor- 
rection to the ground state energy per site SEo /N found in quantum 
Monte Carlo simulations of the quantum ice model [Eq. l|16H. 
Results are shown for cubic clusters of TV = 432, 1024 and 2000 
sites, for parameters yi/g = 0, 0.25, 0.5, 0.75, as function of the 
linear dimension of the system L — {ao/2){N/2)3 . The fact that 
5Eo/N ~ l/L* implies the existence of a linearly dispersing exci- 
tation — the "photon" of the underlying lattice gauge theory. 



In particular, for fi ~ 0, the physical point of our model, we 
find 



(0.6 ±0.1) 5 ao h- 



(116) 



where we have restored the dimensional factor of h. We have 
also calculated an upper bound on c from a single mode ap- 
proximation, in the spirit of Ref. We find 



c < (0.6 ±0.1) gao h'^ 



(117) 



Within errors, the two numbers are indistinguishable. 

It is interesting to use this result to make an order of mag- 
nitude estimate of the speed of light in a quantum spin ice 
material. Considering Yb2Ti207, as (presently) the best- 
characterised material, and inserting the exchange parameters 
obtained by Ross et ai— into the expression for the tunnelling 
matrix element [Eq. (fTsTll. we obtain 



0.05 meV. 



(118) 



From Eq. (11161 ). and the known size of the unit cell 
flo ~ 10.026 A i26ll . we find a speed of light 



c-0.3meVA-50 ms"^ 
which implies a photon bandwidth 

Au - 0.1 meV, 



(119) 



(120) 



within the range accessible to modern inelastic neutron scat- 



tering expenment&sii. 

The accuracy of this estimate is limited by the approxima- 
tions made in setting up the minimal model of a quantum spin 
ice "Htunneiiing [Eq. (fT2]|]. and so it should only be regarded 



as a "ballpark" figure. It should also be remembered that 
Yb2Ti207 is believed to order ferromagnetically at the lowest 



temperature; 



.29.70 



However as long as a given system remains 



an "ice", the inclusion of further tunnelling processes beyond 
^tunnelling should Only increase the speed of light. 



IV. 



'ELECTROMAGNETISM" IN A QUANTUM SPIN ICE 
AT FINITE TEMPERATURE 



In Section UlIBI we have demonstrated that the field the- 
ory Hyjj^-) [Eq. ( l40] il — quantum electromagnetism on a 
pyrochlore lattice — gives an excellent account of the re- 
sults of zero-temperature quantum Monte Carlo simulations 
of the minimal microscopic model of a quantum spin ice, 
Hfj. [Eq. (fTSIll. These results confirm the conjecture that this 
model could support a spin-liquid phase, down to T = 0. En- 
couraged by this, we now use the same field theory to explore 
how correlations in this spin liquid state develop at finite tem- 
perature. 

In Section IIV AI we assess how the thermal excitation of 
magnetic photons changes the temperature dependence of the 
energy-integrated structure factors measured in quasi-elastic 
scattering. We find that pinch-points eliminated by quantum 
fluctuations at zero temperature, are progressively restored as 
the temperature of the spin liquid is raised. 

In Section ITV Bl we compare the results of the lattice field 
theory with published results for quantum Monte Carlo sim- 
ulations of quantum charge ice at finite temperature^. We 
find that both the form and the temperature dependence of the 
correlations are well described by the lattice field theory. 

Finally, in Section IIV CI we conclude with a few remarks 
about the finite temperature behaviour of the heat capacity and 
uniform magnetic susceptibility in a quantum spin ice. 

Throughout this analysis we set ft = fc^ = 1, restoring di- 
mensional factors of h and fc^ only where we quote results for 
the coefficient of the heat capacity associated with photons. 



A. Temperature dependence of structure factors 

The qualitative changes in correlations between spins at fi- 
nite temperature can most easily be understood within the con- 
tinuum field theory Seft [Eq. (l93]ll. The thermal excitation 
of photons enhances correlations of the magnetic field B at 
small Ik I 



X coth 



ck\ 1 



2T 



(121) 



and introduces a thermal de Broglie wavelength for the pho- 
tons. 



At 



(122) 
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FIG. 16: (Color online). Slow, cold death of pinch points in a quantum ice. Equal-time structure factor S^pl„(k,t = 0) [Eq. ( I90H in the 
spin-flip channel measured by Fennell et al\ calculated from the lattice field theory Hyjij [Eq. (I40H. for comparison with neutron scattering 
experiments on a quantum spin ice. Results are plotted for temperatures ranging from T = QioT — ca^^ , where c is the speed of light, and 
ao the linear dimension of the cubic unit cell, with temperature measured in units such that h — ks ~ 1- The pinch-point structure observed 
at finite temperature is progressively "hollowed out" as the system is cooled towards its zero-temperature ground state. 



Over sufficiently long distances, this enhancement of corre- 
lations exactly cancels their suppression by quantum fluctua- 
tions. Assuming that Ac <C Aj, and expanding Eq. (11211) for 
small wave number, we find 

C7,^(|k| « 2n/Xr) = (^6,. - ^) + . . .(123) 

This implies that, for these small wave vectors, the pinch point 
is restored, but with a prefactor that depends linearly on tem- 
perature. 

This result has very simple interpretation. At finite tem- 
perature photons are only coherent quantum excitations over 
a length scale Aj. Therefore, while correlations in a quantum 
spin ice may decay as over distances Ac ^ r ^ Aj, at 
long distances for r ^ Aj the classical 1 /r'' decay of the spin 
correlations is restored. 

All of these arguments generalise to the lattice field the- 
ory "Hfj^-^j [Eq. (l40bl. and to expressions for the equal-time 

structure factor at finite temperatures derived from S"^^^ (k, lu) 
[Eq. (|89tl. Thus we anticipate that they will apply equally 
to a quantum spin ice at finite temperatures. This suggests a 
simple diagnostic for a quantum spin ice in quasi-elastic neu- 




FIG. 17: (Color online) Angle-integrated scattering intensity I{k ~ 
0,r) [Eq. ( I124H calculated from the lattice field theory H'u(i) 
[Eq. ||40H. for comparison with neutron scattering experiments on 
a powder sample of a quantum spin ice. Results are plotted for tem- 
peratures ranging from T = to T = l.Oca^^, where c is the speed 
of light, and ao the linear dimension of the cubic unit cell, with tem- 
perature measured in units such that = fc_g = 1. The progressive 
elimination of pinch points as the sample is cooled manifests itself 
as a steady loss of scattering for |k| — > 0. 



tron scattering experiments — as the sample is cooled, and 
photons become coherent over longer length scales, the pinch 



22 



points observed at reciprocal lattice vectors are progressively 
"bleached out". This slow, cold, death of pinch points is illus- 
trated in Fig.[T6l 

Since there is also a characteristic loss of spectral weight 
in >5'^fn(k, t = 0) for k w 0, exactly the same process could 
be seen in the angle integrated structure factor measured in 
neutron scattering experiments on powder samples. In this 
case, the intensity of scattering is given by 



a/3 ^ 



(124) 



For classical spin ice, or a quantum spin ice at sufficiently high 
temperature, 

I{k « 0, T) « const. 

However, as a quantum spin ice is cooled to zero tempera- 
ture, the growing coherence of photons will manifest itself as 
a progressive loss of spectral weight at small k, 

I{k = 0,T) - T 

until, for T = 0, spectral weight at fc = is eliminated en- 
tirely 

/(fc w o,r = 0) cx fc 

This progression is illustrated in Fig. (TT] 

B. Comparison with quantum Monte Carlo simulation 

It is also interesting to compare the predictions of the lat- 
tice field theory H'^j^^i^ [Eq. (l40bl. with the results of finite- 
temperature quantum Monte Carlo simulations of a quantum 
charge ice described by "Ht-v [Eq. (flTtl. as published by 
Banerjee et al^. Banerjee et al. performed their simulations 
for hard-core bosons on a pyrochlore lattice at half filling, 
with hopping integral t = 1, and nearest-neighbour repulsion 
V = 19.4, at temperatures T = l.OSg and T ~ 1.57g, where 
g = 12t^ /V^ is the size of the leading tunnelling matrix ele- 
ment between different charge ice configurations. 

In Fig. [T8]we plot simulation results for S'charge(k, t = 0) 
at these temperatures, calculated within a single sublattice of 
pyrochlore lattice sites, together with the best fit to Eq. ( 1831 ). 
projected onto a single sublattice. We assume that the param- 
eters of the field theory depend relatively weakly on tempera- 
ture, and attribute the temperature dependence of correlations 
entirely to the thermal excitations of photons. Under these as- 
sumptions, the lattice field theory gives a good account of both 
the form and the temperature dependence of Schargelk, t = 0), 
within the error bars on points taken from simulation. 

These fits suggest a speed of light 



(1.8 ±0.2) gao fT 



(125) 



which is '-^ 3 times larger than that found in Section lTlI DI from 
finite size scaling of the ground state energy of T-L^ [Eq. ([T6bl. 



This discrepancy can probably be attributed to the fact that 
the simulations of Banerjee et al. were performed close to the 
melting point of the charge ice^, where both interactions be- 
tween photons, and tunnelling processes involving more than 
six lattice sites, are likely to play an important role. Since 
all of these processes will contribute to the rate at which the 
gauge field fluctuates in time, they can be expected to increase 
the speed of light. 
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FIG. 18: (Color online). Comparison of the predictions of the lattice 
field theory H'u(i) [Eq. ( I40H with the results of finite-temperature 
quantum Monte Carlo simulations of a quantum charge ice described 
by Ht^M [Eq. J17H. Results are shown for the equal time, on- 
sublattice structure factor 500(14) = (no(— k)no(k)) Simulations 
are taken from Banerjee et al.^^, and were performed at tempera- 
tures r = 1.05^ and T = 1.57g, where g = 12*^/1^^ is the size 
of the leading tunnelling matrix element between different charge ice 
configurations. The temperature dependence of the spin correlations 
makes it possible to estimate the speed of light c ~ 1.8 g ao h"^ . 



C. Heat capacity and magnetic susceptibility at low 
temperatures 

Neutron scattering experiments have the potential to give 
decisive information about emergent electromagnetism in a 
quantum spin ice. However these experiments are expen- 
sive and difficult to perform, and depend critically on the 
size and quality of available samples. We therefore conclude 
with a few brief remarks on potential signatures of a quantum 
U{1) liquid in thermodynamic quantities. The results given 
will hold in the low-temperature regime where the physics 
of a quantum spin-ice can be described as a gas of photons. 
At higher temperatures the thermal excitation of the gapped 
spinons (monopoles) and electric charges also play an impor- 
tant role. 

We have seen in Section IIV Al how quantum fluctuations 
lead to an equal-time structure factor which, in the limit 
k ^ Q, vanishes at low temperatures as 



lim S'(A:, T) (x T 



(126) 
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This in turn implies a bulk magnetic susceptibility x(T) which 
is independent of temperature at low temperatures 



^-\T«g) = ^. 



(127) 



where U is the coefficient of in the effective Hamiltonian 
^u(i) [Eq- (l38])l, and k w 1 is the dimensionless scale factor 
introduced in Eq. ( f74b . 

This result provides another means of parameterizing the 
lattice field theory. It is also a potentially useful diagnostic 
for experiment, since, a classical spin ice which remains in 
thermodynamic equilibrium at low temperatures, will exhibit 
an effective Curie lav 



(128) 



11.81 



This result follow directly from the fact that there are more 
spin ice configurations with vanishing magnetisation M = 
than with any finite magnetisation M 7^ 0, and so, in the ab- 
sence of any other considerations, a state with M = is se- 
lected by an entropic term 5F ~ T5S ^ T in the free 
energy^. Nonetheless, any comparison with a classical spin 
ice should be approached with some caution, as these sys- 
tems need not remain in equilibrium at low temperature s ^^'^-^ , 
and the character of the spin fluctuations which control x(r) 
changes as a function of temperatur e 

As noted elsewhere^ilS, the fact that photons are linearly 
dispersing excitations implies that they must make a con- 
tribution to the heat capacity at low temperatures. While this 
contribution has exactly the same temperature dependence as 
that from acoustic phonons, the large amount of entropy avail- 
able in ice states, and low speed of light [cf. Section UlI PL 
mean that the heat capacity at low temperatures will be domi- 
nated by photons. The photon contribution to the heat capacity 
per mole is 



Cphoton[™o/e] = BT^ 



with the coefficient B given by 



B 



R 



he 



(129) 



(130) 



From the characterisation of Yb2Ti207 by Ross et al^, and 
the analysis of the speed of light in Section ITlIDI we estimate 



B 



65 J mol"^ 



(131) 



which is several orders of magnitude larger than the expected 
phonon contribution. This should be compared with the value 



1 J mor^ K"'' 



obtained in Ref. |70l We note that, since the photons are mag- 
netic excitations, measurements of the heat capacity in an ap- 
plied magnetic field may also prove instructive. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have developed a detailed theory for 
the simplest microscopic model which could describe quan- 
tum tunnelling between different spin ice configurations 
[Eq. (fT2bl. The striking claim that this type of model 
could support a spin liquid phase which perfectly mimics 
quantum electromagnetismii has been verified by quantum 
Monte Carlo simulations^^iS. Here we have explored how 
such a quantum spin liquid might manifest itself in experi- 
ment, parameterizing an "electromagnetic" lattice gauge the- 
ory from quantum Monte Carlo simulations at zero tempera- 
ture, and using this to calculate the dynamical structure factor 
5"'^(k, w) [Eq. ( |89] l1 which would be measured in neutron 
scattering experiments at finite temperature. 

We find that a key signature of the emergent electro- 
magnetism is the suppression of pinch points singularities 
in the energy-integrated structure factor S'"''(k, < = 0) as 
the system is cooled to its zero-temperature ground state 
[Fig. [T6l . This will coincide with the appearance of a gap- 
less, linearly dispersing, mode — the photon of the lattice 
gauge theory — in inelastic neutron scattering [Fig. [T3l . 
In sharp contrast with a conventional antiferromagnet, the 
dispersing feature associated with this photon vanishes as 
uj ^ 0. These photons will also strongly influence the low- 
temperature thermodynamic properties of the system, giving 
rise to a temperature-independent contribution to the magnetic 
susceptibility [Eq. (11271 )1 and an anomalously large contri- 
bution to specific heat [Eq. ( |129H ^. 

Neither the idea of "artificial light"— nor the observa- 
tion that there could be quantum tunnelling between different 
spin ice configurations^"", is new. However the possibility that 
one might lead to another is both new and exciting, and adds 
to the general frisson surrounding pyrochlore magnets. With- 
out attempting to review all of this fast-developing field — but 
with the possibility of observing photons in mind — it is in- 
teresting to ask whether any of the materials currently studied 
"fit the bill". 

The most widely studied example of a three-dimensional 
spin liquid is the insulating pyrochlore oxide Tb2Ti207^. 
Tb2Ti207 does not order down to 50 mK 1231] . despite hav- 
ing a Curie- Weiss temperature 6cvj ~ 14 K jlSII . and a strong 
tendency to order under magnetic field or pressur e In 
a series of papers, Gingras and coauthors have argued that 
Tb2Ti207 is a "quantum spin ice", in which spins fluctuate 
strongly about the crystallographic [111] axes. These claims 
were made on the basis of a characteristic checkerboard struc- 
ture observed in diffuse neutron scattering experiments at high 



temperature! 



,19.88 



and a subsequent microscopic analysis of 



crystal field levels^Si^, and find support in the recent observa- 
tion of partial magnetisation plateau for magnetic field applied 
along a [111] axis2ii22. 

Within this framework, the field at which the plateau is ob- 
served implies that the energy scale relevant for "quantum 
spin ice" behaviour in Tb2Ti207 is Jeff « 0.2A=. Unfor- 
tunately, the interpretation of experiment at these low tem- 
peratures is muddied by questions of sample quality, with in- 
consistent results for spin-freezing obtained by different au- 
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thors22i212ii22. Published thermodynamic data at low temper- 
atures is also less than conclusive, showing hints of a satu- 
ration of x(T) at low temperatures, but strong sample de- 
pendence2i22i^. And the picture is further complicated by 
strong fluctuations of the lattic o^^'^^ , with alternative theories 
of TbsTizOr building on lattice effects^Sriffi 

At present there is no published neutron scattering data for 
Tb2Ti207 with the combination of k-resolution, energy reso- 
lution and low temperature needed to compare with the pre- 
dictions in Section HTD] and Section ITlI CI of this paper How- 
ever recent evidence of "pinch-point" structure in quasi-elastic 
neutron scattering on single crystals of Tb2Ti2072^, taken to- 
gether with inelastic neutron scattering experiments on pow- 
der samples^, suggest that the comparison might be interest- 
ing. The latter find evidence of a quasi-elastic feature evolving 
into two bands of excitations at temperatures T < 0.4Kii2ii. 
If — and it remains a big IF — the behaviour of Tb2Ti207 
is connected with the physics of the quantum ice described 
in this paper, it would be tempting to identify these bands 
with the excitations of electromagnetism on a lattice — gap- 
less photons, together with gapped "electric" and "magnetic" 
charges (spinons). But more, and more delicate, experiments 
will be needed to determine whether this is indeed the case. 
And ultimately, Tb2Ti207 will remain a fascinating system to 
study, regardless of whether or not it is a quantum spin ice. 

Recently, there has also been intense experimental and 
theoretical interest in the closely -related Yb pyrochlore, 
Yb2Ti207. Originally identified in the pioneering survey 
of Blote et alJ^ as a ferromagnet with Tc = 0.21K and 
dcvj = 0.4K, Yb2Ti207 differs from the classical spin ice 
materials Ho2Ti207 and Dy2Ti207 in that the lowest ly- 
ing crystal field state is a Kramers doublet with easy-plane 
anisotropy^2iSiS. An XY ferromagnet on a pyrochlore lat- 
tice — m odern estimates suggest 9c\n ~ 0.65Kfor Yb2Ti207 
ll67lll06ll — would naturally be expected to order ferromagnet- 
ically at low temperatures. However Yb2Ti207 exhibits a far 
more complicated phenomenology. 

Neutron scattering experiments at temperatures below lOK 
find diffuse liquid-like structure which offers evidence of 
anisotopic exchange interactions^^^. At temperature of order 
IK, rod like structure emerges, reminiscent of a dimensional 
crossove r~^'^^' ' "^- .los ^ Some authors have found evidence of a 
first order transition into a ferromagnetic ally ordered state at 
Tc = 0.24K22ii2^, although this has been contested, and may 
not occur in all samplea''^-'°^i^"^i'°A"'s. 

That Yb2Ti207 orders ferromagnetically in applied mag- 
netic field is, however, uncontroversial. And this has made 
it possible for Ross et alM- to accurately characterise an ex- 
change Hamiltonian for Yb2Ti207 [Eq. ( flOt l from fits to spin 
wave excitations in the polarised state. The parameters ob- 
tained confirm that the dominant interactions in Yb2Ti207 
favour "ice" states, but that these are complimented by terms 



which will drive significant fluctuations at low temperatures. 

Reassuringly, this description of Yb2Ti207 is also in 
quantitative agreement with measurements of thermodynamic 
properties over a wide range of temperatures^. This makes 
Yb2Ti207 the best-characterised "quantum spin ice", and as 
such, it is a natural place to look for emergent electromag- 
netism. However neutron scattering data with sufficient reso- 
lution to compare with the predictions of this paper are not, as 
yet, available. 

Tb2Ti207 and Yb2Ti207 are by no means the only py- 
rochlore systems with spin-liquid properties^, and some 
of these other systems, notably Pr2Sn20 | ^'-^^i^o^ and 
Pr2Zr207>^^ are also worth investigating as potential real- 
isations of a quantum ice. It might also be interesting to 
revisit two-dimensional ice-type materials, such as the pro- 
ton bonded ferroelectric copper formate tetrahydrateii. While 
two-dimensional quantum ice models are known to order at 
low temperature sJJ^—, they are described by the same class 
of lattice gauge theory, and possess the same spinon excita- 
tions as their three-dimensional counterparts^SiliLli^. These 
excitations will be confined in the ordered state, but might be 
visible at finite energy, and above the ordering temperature . 

Although the theoretical possibility of emergent elec- 
tromagnetism in quantum icelil2i^2i5S and quantum 
jjjjjgj46j60j6ij72 models is now well-established, many 
theoretical questions remain open. In this paper we have 
considered only the simplest microscopic model of a quantum 
spin ice [Eq. (fT6bl, and fully characterised only its photon 
excitations. The study of more realistic models, and of other 
excitations, is still in its infancy^2i22ili. We have also made no 
attempt to resolve the question of how the quantum ice state 
which we find at low temperatures, becomes a classical ice 
at high temperatures. All of these issues remain for future 
study. But we believe that the best motivation for studying 
them is experiment, and hope that the results in this paper 
will encourage further experiments on spin liquid materials 
which might realise artificial light. 
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